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Abstract. We extend a linear version of the liquidity risk model of ^etin et al. (2004) to allow 
for price impacts. We show that the impact of a market order on prices depends on the size of the 
transaction and the level of liquidity. We obtain a simple characterization of self-financing trading 
strategies and a sufficient condition for no arbitrage. We consider a stochastic volatility model in 
which the volatility is partly correlated with the liquidity process and show that, with the use of vari- 
ance swaps, contingent claims whose payoffs depend on the value of the asset can be approximately 
replicated in this setting. The replicating costs of such payoffs are obtained from the solutions of 
BSDEs with quadratic growth and analytical properties of these solutions are investigated. 

1. Introduction 

In financial markets, liquidity either refers to the ease with which financial securities can be 
bought and sold or to the ability to trade without triggering important changes in asset prices. 
Liquidity becomes a risk factor when the magnitude of the impact of these phenomena changes 
randomly over time. Uncertainty regarding the level of liquidity in traded assets has been for a 
long time a critical issue for moderate to large traders. The cost of a given trading strategy in real 
world situations can be substantially high when large quantities of financial assets are traded due 
to the consequential impact of trading on prices, and the limited and uncertain future supply and 
demand. In this paper, we construct an arbitrage-free model which relates levels of liquidity to 
trade impacts and quantify liquidity costs of strategies used for hedging claims contingent on the 
value of the traded asset. 

The literature on liquidity risk is large and can be mainly divided according to these two con- 
ceptual perspectives. In the first category of models, the price of an asset depends on the size of the 
transaction and the depth of the order book. The second category includes those commonly known 
as "large trader" models in which a large trader buys and sells such large quantities of assets that 
his trades affect the prices in a non-negligible way. The purpose of this paper is to combine both 
approaches in a unified framework and to study the problem of contingent claim replication. 

Examples of recent papers in the first category include ^etin, Jarrow and Protter [|6l and ^etin 
and Rogers Q. Rogers and Singh [|20l give a microeconomic argument for a price which depends 
on size and this is then reflected in the dynamics of self-financing strategies. They solve an optimal 
control problem in this context. 

Bank and Baum BH, Frey [12J and Jarrow [fTSl are examples of papers in which the impact of 
the large trader is a function of its current holdings. In Alfonsi et al. [HI, the authors relate the 
impact of trades to the shape of the order book and consider the problem of optimal liquidation 
by the large trader. On the other hand, Ly Vath et al. [17J study the problem of optimal portfolio 
selection for a large trader who has a price impact function and cost function of exponential form. 

Our present model was in part inspired by the liquidity risk model of ^etin, Jarrow and Protter 
(61 (thereafter referred to as the CJP model). In the CJP model, liquidity is introduced by hypoth- 
esizing the existence of a supply curve S{t,x) which gives, at a given time t, the price per share to 
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pay for a stock in terms of the size x of the trade. In such a model, the trader observes the supply 
curve and acts as a price taker. In this setting, liquidity costs essentially depend on the quadratic 
variation of the trading strategy. The main drawback of this model is that liquidity risk can essen- 
tially be avoided by approximating a given self-financing trading strategy (s.f.t.s.) by a sequence 
(Z")„>i of continuous s.f.t.s. with finite variation (FV) which incur no liquidity costs. The prices 
of options are then unaffected by liquidity risk. This issue was cleverly dealt with in ^etin et al. 
[[8 1 by adding constraints on the gamma of the hedging strategies. A liquidity premium is then 
reflected in option prices. 

Our approach is to combine both notions of liquidity risk by hypothesizing the existence of 
a random linear supply curve and by studying the impact of trades on prices. One of the key 
observations made in this paper is that the magnitude of price impacts is directly related to the 
amount of liquidity of the asset. This leads to a simple characterization of self-financing trading 
strategies in which the profit is partly affected by the level of liquidity. The main goal of this paper 
is to study the effect of liquidity risk on the replicating costs of contingent claims. We consider a 
stochastic volatility model in which the volatility process depends in part on the level of liquidity. 
We will see that variance swaps are the simplest hedging tools in this setting. 

The paper is organized as follows. In Section |2l we derive the impact of trading on prices using 
simple principles and show that changes in the price of an asset is directly affected by the amount 
of liquidity. We then use these observations to propose a model defined on the Brownian filtration 
and show it is arbitrage-free. A simple characterization of self-financing strategies is derived to 
help set up the replication problem. Section [3] is devoted to the main result of this paper, the 
replication of contingent claims using variance swaps and the characterization of replication costs 
in terms of backward stochastic differential equations with quadratic growth. Section |4] presents 
useful analytical properties of these solutions. 

2. The Setup 

We consider an economy consisting of a risky asset (typically a stock) which is traded through 
a limit order book, its associated contingent claims and a risk-free asset. We take the point of view 
of a hedger who observes the limit order book of the stock and makes market orders (also known 
as marketable limit orders). We start by describing the supply curve the hedger would expect to 
observe if he did not trade. We call it the unaffected supply curve and denote it by S. It represents 
the limit order book that results from all other traders' limit and market orders. It is a conceptual 
construction which is not directly observed. We will assume that the hedger's trades have a lasting 
impact on prices which will be added to S to obtain the actual observed supply curve, which we 
denote by S'^. 

We are given a fixed maturity T and (Q, ^ , (^?)o<?<r5 P) a filtered probability space which sat- 
isfies the usual conditions. We assume the interest rate is constant, and for simplicity we work with 
the discounted price processes. The (discounted) unaffected price process is an exogenously given 
adapted continuous process 5= (5'f(x))^>o,xeK (or sometimes written 5 (?,jc) for convenience). St{x) 
is the price per share for a transaction of size x at time t that would be observed if the hedger did 
not trade before time t. The actual (discounted) quoted price per share that all market participants 
obtain for a transaction of size x at time t is denoted by S^{x). We start by assuming that the 
unaffected supply curve has the following linear structure: 



Sf{x) = St +Mtx, for xeR 
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(2.1) 



where {Mt)t>o and (St)t>o are positive semimartingales. Note that the fact that this function is 
continuous atx = implies there is no bid-ask spread. While it is theoretically possible for St{x) to 
be negative for some values of x, it is unlikely to happen in practice since the value of Mt is small. 
We assume there is a measure Q, equivalent to P, under which the unaffected price process S is 
local martingale. As in the classical theory, this assumption will be sufficient to rule out arbitrage 
opportunities. See Theorem 12 . 5 1 below in this regard. 

The assumption that the supply curve is linear is supported by the empirical study of Blais f3l 
for stocks that are frequently traded in large volumes. The study was based on a large data set of 
stocks traded on several different stock exchanges in the year 2003. See also Blais and Protter [4J. 

Before we specify the precise model for S and M on which we will focus, we start by detailing 
general characteristics that a liquidity risk model which include price impacts should reflect. 

Equation (12.11 ) gives us a way to describe the limit order book. We represent it by a density 
function pt{z) which denotes the density of the number of shares being offered at price z at time t, 
i.e. J^^ pt{z)dz is the number of shares offered between prices zi and zi- If a trader wants to buy 
X shares at time t through a market order then the price to pay is J^" zpt{z)dz in which Zx solves 
Ss^ Pt {z)dz = X. It is clear from the linear structure of the supply curve that for any t <T the density 
equals pt{z) = In that case, Zx = 5"? + IMfX and the dollar outlay for x shares is 



Since p is a measure of liquidity, we can think of M as a measure of illiquidity. Indeed, the larger 
is Mj, the higher is the liquidity cost. 

We let Xt denote the number of shares owned by the hedger at time t and S^{x) denote the actual 
asset price per share observed in the market, which includes the impact of the hedger's trading 
strategy, i.e. ^{x) implicitly depends on X. We define 5f = S^{0) as the observed quoted price. 

We now describe the impact that an arbitrary market order has on the limit order book. We will 
then use these observations to justify our specification of S and S^. First, one should observe that 
if hXt shares are bought at time ? by a trader through a market order, then the corresponding part 
of the order book is used up. This would mean that immediately after the trade the limit order 
book would have a density of for prices between and 5f + 2M/AX, and pt elsewhere since the 
lowest ask price would then be + 2M/AX, whereas the highest bid would remain the same. In 
this perspective, one can see that an implicit assumption made in the liquidity model of ^etin et al. 
|l6l is that new limit orders to sell are placed immediately after a trade, thereby filling up the limit 
order book to its previous levels since it is assumed that trades have no impact on the supply curve. 
The new observed quoted price is the same as before and the impact on prices is non-existent in 
li6l . Although it is reasonable to assume that the limit order book fills up to its previous level after 
a trade, it is not clear whether the gap should be filled by bid or ask orders. For example, if the 
gap is filled entirely by bid orders after the purchase of AXf, then the new quoted price is shifted 
upwards to + 2A//AX/. In this case, the outcome is a full impact on prices. 

The empirical findings of Weber and Rosenow [2T| showed that in practice the impact of trading 
on prices is important but can be less than the full impact described in the previous paragraph. In 
fact, they showed a negative correlation between returns and the volume of incoming limit orders 
which suggests that traders respond to buying market orders by adding new limit orders in the 
opposite direction. We model this phenomenon by introducing a parameter A G [0, 1] measuring 
the proportion of new bid orders (resp. ask orders) filling up the limit order book when a trade to 
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buy (resp. sell) is made at time t. In effect, the effective impact on prices of a trade of size AX/ is to 
shift the quoted price to 5f + 2AM/AX/, whereas the density level of the order book is unaffected. 

We have to be careful how we define the observed price process in this setting. Indeed, when 
the hedger makes a trade at time t the price he pays is unaffected by the impact of this current trade 
whereas prices right after t will be. In this sense, will not be cadlag in general, although S^j^ is 
and includes the impact of a trade at time t. 

Suppose X is a simple trading strategy of the form X = X^^'^q A^X1[^«^-) in which A^X = X^-^" — 
X^n ^ for k= 1 , . . . , n and AJjX = Xq. Then, the observed quoted price should satisfy 

k-\ k-l k-l 

S^ = St + 2Y, ^M^nt^lX = Sf + 2 £ XM^n_^l^]X + 2 £ A (AfM) (A^X) 

for any t E 'T'^']. Note that the sum in the previous equation only goes up to fc — 1 since 5°„, 

which represents the price per share for a trade of size 0, is not yet impacted by the trade at time 
T'j^. The right-limit version of this process is then given by 

k-l k-l 

50^ = 5^ + 2 £ AM^«_ Af X + 2 £ A ( Af M) (Af X) (2.2) 

;=0 ' 1=0 

for any t G 'T'^) when S is right-continuous. Following these observations, we define 

5«+ = S, + 2A f Mu-dXu + lX f d[M,X]u (2.3) 
Jo Jq 

for all t <T, for a general semimartingaleX. Furthermore, we define the observed quoted price by 
= lim^if 5^_|_ . By assuming that the level of liquidity pt is unaffected by trades, we readily obtain 
that the supply curve is given by 

S^{x)=S^ + M,x (2.4) 

for all < f < r and jc G M. We think of 1 — A as the fraction of the order book which is renewed 
after a market order so that in practice the actual impact on prices is A times the full impact. 

Equation l2.3l gives us a new understanding of the causes of volatility and its relation to illiquidity. 
As mentioned earlier, S is the price process which results from limit and market orders of all the 
other market participants. The equation suggests that the impact of the market orders of each 
market participant is proportional to the value of M. The volatility of S can then be expected to be 
correlated in part to M. (Another component of the volatility of S would be related to the volatility 
of limit orders.) In fact, many empirical works have shown that the level of liquidity is an important 
determinant of the variance of log-returns. The reader is referred to the works of Farmer et al. [11 J 
and Weber and Rosenow [|22l for a more detailed discussion. The observation that these authors 
make is that volatility is high when liquidity is low, and low when liquidity is high. Since M is 
a measure of illiquidity, we can expect the instantaneous variance of the log-returns of the stock 
price to be in part correlated with M. This is a key observation which will enable us to hedge 
derivatives. Indeed, in the next section, we will introduce variance swaps which will be used to 
hedge against the liquidity risk. Since volatility is one of the most correlated quantities to liquidity 
risk, this is a very natural choice. See Remark |3] in this regard. 

Following these observations, we consider a stochastic volatility model for S: 



dSt = 'LfSfdWij, 
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(2.5) 



in which Wi is a Brownian motion defined on the filtered probability space, and Ef is the stochastic 
volatility. Recall that we are working directly under a risk neutral measure Q for unaffected prices, 
hence S has no drift term. We model M and E as follows. Define V and U as the solutions of 

dUt = riUt + ri)dt + ^{Ut)dW2,t, 

dVt = a{Vt + a)dt + &{Vt)dW3^t 

in which W = (W^7,?);<3,f<r is a three dimensional Brownian motion defined on the filtered proba- 
bility space, and a, 7,7], a G M. We define = Ut + Vt and let M = £r{U), in which F is strictly 
increasing and twice continuously differentiable. In practice, the process M takes small values 
compared to E, but is also an important component of the volatility process E. As a result, the 
constant e is typically small. 

We are using a three dimensional Brownian motion since there are three different sources of risk 
in this model, namely the stock price, the liquidity level and the volatility, which is, in practice, 
only partially dependent on the level of liquidity. The components of W are typically correlated 
and we denote by R = jCOV{Wt) the matrix of instantaneous correlation coefficients. We assume 
R is positive definite and we let L be the upper triangular matrix in the Cholesky decomposition 
such that = LJL. We then define B = LW. Then 5 is a three-dimensional Brownian motion 
with independent components. We denote the components of by 

h h \■ 
03 J 

We assume the functions and 4> are chosen so that the solutions of the above stochastic dif- 
ferential equations are well defined. For example, one can take 0(v) = with = 0,^ or 1. 
Examples of stochastic volatility models of this form are Heston [13l (6 — Hull and White [14J 
(0 = 1), and Detemple and Osakwe fW\. Other expressions for I? could be used, however we 
have chosen this particular form for its mathematical tractability and its widespread use in theory 
and practice. 

2.1. Self- Financing Strategies and No Arbitrage. In order to properly address the problem of 
replicating contingent claims, we give a characterization of self-financing strategies and establish 
under which condition our model is arbitrage-free. In our setting, the self-financing condition is as 
follows. 

Definition 2.1. Let 7t„: to = Tq < r" < . . . < T^'^ = T be a sequence of random partitions tending to 
the identity. A pair of processes {Xt,Yt)tQ<t<T is a self-financing trading strategy (s.f.t.s.) on [to, T] 
if X is a semimartingale and Y is an optional process satisfying 

Yr = y,„_-AX,o50(?o,AX,J-limX;A^X50(T,",A^X)l|,„<,} (2.6) 

where convergence is in ucp. (See Protter |fT9ll for undefined terms.) Here, A'^X = X^n —X^n for 
k= 1, . . . ,n. 

Xt represents the number of shares of the asset owned by the hedger and Yf is the position in the 
risk-free asset at time t. The interpretation is that the position in the risk-free asset at time t should 
be equal to the position at time to minus the cost of all the trades between to and t. Here, 7/^- is the 
value of the position in the risk-free asset before the trade at time Iq. 

5 



Remark. In the classical theory, the process X is predictable. We take X in the above definition to 
be a semimartingale for the stochastic integral in Equation 12. 3 1 to be well defined. A consequence 
of Proposition 12.21 below is that the limit in Equation 12.61 is well-defined, and the definition of 
self-financing trading strategies is independent of the sequence of random partitions. 

Even though s.f.t.s. are defined in terms of 5°, they can be characterized in terms of the exoge- 
nously given processes M and S as follows: 

Proposition 2.2. Let to > 0. If {Xs,Ys)tQ<s<T is a self-financing trading strategy then 
Yt+XtiS'i^-XMtXt) = Yt,-+Xt,-{Sl-XMt,Xt,-)+ fxu-dSu 

Jto 

-A fxl_dMu - [\l -A)M„J[X,X], (2.7) 

Jto Jto 



for all to <t <T. 

Proof. Let Ttn : to = Tq < < . . . < x'^ = ? be a sequence of random partitions tending to the 
identity. The self-financing condition is 

k 

Y, = y,o--^fo^°(^o,AX,J-lim£Afxte+M,«Afx) 

i=\ 

where the convergence is in ucp. We can expand the sum in the last equation to find 

- iim„^<.i:f:Li ^"^ {fx'i +^T«Afx^ 



= - lim y fc5°„ -X^n ) + lim y X^n A"S^ - lim f MA^'X 

1=1 1=1 1=1 

= + +Jim (X,«_A'^S + 2XM,«_X,>._AUX^ - Jim £M,«(AfX^2 

;= 1 i= 1 

k 

= -Xf 5° + Xr.sl + 2XMt,Xt,AXt, - iXMtXfAXt + lim f X,« , A;^5 

!=1 

k k k 

+2 lim y\XM^nX^nAlX- lim V^XMA^lXf- lim V (1 - A)M^n(AfX)2 

;= 1 != 1 (= 1 

k k 

= + XtX+ + lim £ X,n Af 5 + lim f XM^nA^X^ 

i= 1 != 1 

k 

-lim£(l-A)M,«(AfX)2 

!=1 



= -XtS^+ + Xt,Sl+ + XM,X^ - XMt,Xl + lim £ X,n Af 5 



lim y XX^n I^IM 

Yl — >oo /— 1 




i=l 



-limy(l-A)M^«(A"X)2 



7=1 




-A / Xl_dMu - / (1-A)M„J[X,Z] 



by Theorem 21 (Chapter II) of Protter lfT9ll since X is cadlag. 



□ 



One can think of Yt+x{Sf — XMfx) as the liquidation value of a portfolio with x shares at time t. 
Indeed, take tQ = t and AXf = X/ in Equation l2.7[ Then one finds that the cash value of a position 
Xt- at time t- in the stock is equal to AF, = Xt- {S^ - XMfXt-) - ( 1 - A )Mf (AZf )2 if it is liquidated 
at time t. Furthermore, if one uses a sequence X" a continuous and FV processes converging to 
X (this can be done by Lemma [X2l ). then the liquidation value converges to Xt-{S^ — XMfXt^). 
Consequently, XMt can be interpreted as the effective liquidity parameter. 

Similar to the infinitely-liquid case (M = 0), Equation 12.71 states that the difference in the liq- 
uidation values between time Iq and t is equal the cumulative gains in the risky asset Jl^^Xu-dSu, 
except that in this case there are added costs coming from the finite liquidity of the asset. First note 
that if A = we get a linear version of the CJP model. The integral with respect to M is related 
to the impact of trading. If A = 0, the limit order book is automatically refilled after a market 
order, as in the CJP model. At the other extreme, when A = 1 the impact of trading is at its fullest. 
It is interesting to notice that whatever the trading strategy used an investor always has a partial 
benefit from the asset becoming more liquid. Indeed, when Mt decreases, the associated integral is 
positive no matter what the sign of Xt is. To understand this, it is important to remember that the 
hedger's trades have a permanent impact on the quoted price which is proportional to the level of 
liquidity. If the liquidity is low when he purchases a share and high when he sells it, the price goes 
up higher after his purchase then it comes down after the sale. As a result, the hedger has a partial 
gain from this trade. This is a typical characteristic of large trader models. Note that unless the 
hedger uses a trading strategy with zero quadratic variation this is only a partial benefit because 
there is always a liquidity cost associated to his trades. 

Using Proposition 12.21 for j G M, we define the set of payoffs of maturity T attainable at 
price y by ^7 -measurable random variables Yt of the type 



in which (X,),>o is cadlag with finite quadratic variation. 

def 

We will denote by ^ = UyeK the set of all attainable payoffs. We use the following defini- 
tion of admissibility. 

Definition 2.3. Let a > 0. A s.f.t.s. {Xt,Yt)t>Q is a-admissible if 



y+ Xt-dS,-X Xf_dMt~ {I - X)Mtd[X ,X]t 







Xs-dSs-X / X]-_dMs- / {l-X)Msd[X,X]s>-a 
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for all t <T . The s.f.t.s. {Xt ,Yt) s^t>Q} is simply said to be admissible if it is a-admissible for some 
a>0. 

A strategy is admissible if its payoff is bounded from below. In particular, this definition rules 
out doubling strategies and is well known to be a key element in the definition of arbitrage oppor- 
tunities. See Delbaen and Schachermayer [9] in this regard. 

Definition 2.4. An arbitrage opportunity is an admissible s.f.t.s. whose payoff Yj E satisfies 

F{yr > 0} = 1 and ¥{Yt > 0} > 0. (2.8) 

It is already known (see f6l) that the existence of a local martingale measure for S rules out 
arbitrage opportunities in the CJP model. In the presence of trade impacts, the equation for the 
payoff of a s.f.t.s. has an integral with respect to M. Since the integrand of this integral is always 
negative {—XX^_), then the part of the profit coming from this integral will be negative on average 
if M is a submartingale under the risk neutral measure. This idea is made precise in the following 
theorem which gives a sufficient condition for no arbitrage. 

Theorem 2.5. If there exists a measure Q ~ P under which S is a Q-local martingale and M is a 
Q-local submartingale, then there are no arbitrage opportunities. 

Proof. By the Doob-Meyer decomposition theorem there exists a Q-local martingale M and an 
increasing predictable process A such that M = M +A. Let Zf = Jq Xu^dSu — A /q X^_dMu — /q ( 1 — 
X)M,d[X,Xl. Then Z, + A Jl^X^^dA, + Jl,{l - A)M„rf[X,X]„ = J^X.^dSu - A JoX^^dM^ > -a 
since A is increasing and /o(l — X)Mud[X,X]u > 0. Now, S and M are Q-local martingales hence 
Zf -f A /qX^_JA„ + /o(1 — X)Miid[X,X]u is also a local martingale and because it is bounded from 
below it is a supermartingale. Therefore, Z is also a supermartingale and EqZj- < 0. But, because 
Q ~ P, if Zt were an arbitrage opportunity it would also satisfy Equation 12 81 with Q instead of P 
and EqZt > 0. □ 

In the simplest case, when r{x) = x, it suffices to take y and r] > 0. In the case that r{x) = x^, if 
4>(m) = y/Tn then we need yr] > ^; if 4>(m) = m then we must have 7 > ^. In the remaining parts 
of the paper, all expectations are with respect to Q. 

3. The Replication Problem 

We now turn to the problem of contingent claims replication. Because the presence of the 
processes M and E involve risks that cannot be hedged completely by solely trading the stock, not 
all payoffs are attainable when only the underlying asset is allowed to be traded. Because these two 
processes are components of the instantaneous variance of the log-returns of the stock, the natural 
hedging instruments to consider are variance swaps. We thus consider contingent claims denoted 
by Gi {i = 1 , 2) for which the payoff at time Ti> T {T\ ^ T2) equals the difference between the 
realized variance over the time interval [0, 7]] and a strike Ki, i.e., 

Gij. = l^'i:]ds-K,= C'{Us + Vs)ds-K,. 
Jo Jo 

To rule out arbitrage opportunities, we assume the unaffected price processes G' are Q-martingales 
(z = 1,2), i.e. 



for all t < Ti. We further assume the G/'s have a linear supply curve, i.e. Gi^t{x) — Gi^t for all 

X and t <T . Since it is not infinitely liquid, trading can affect its price and we denote by A/M- ^ 
its effective liquidity. Typically, changes in the supply curves of the G/'s will happen less often. 
Hence, to keep the problem tractable, we assume that M- ^ = M- q is some given positive constant 
for all ? e [0, r] . We will see that two of these swaps are sufficient to hedge against liquidity risks. 
Because we now have two more traded assets, X\^t denotes the number of shares of G\ and X2,t the 
number of shares of G2 in the portfolio at time t. We can easily extend the definition of s.f.t.s. to 
the case of three traded securities. As shown before, s.f.t.s. (X,;^,y) satisfy 

Yt = Yt+ f Xu-dSu+Y. C Xi,u-dGi,u-^ C xl_dMu 

- [\\-X)Mud[X,X]u- £ l\l-K)M[j[XuXi]u 

-^t i=l,2 ' 

for to<t<T, when = XUo = X2,to = = Xi,T = X2,t = 0. 

The following proposition gives condition under which the three price processes 5, Gi , G2 are 
non-redundant. It justifies the choice of variance swaps by providing a simple explicit representa- 
tion of the processes G/. This result will then be used to solve the replication problem. 

Proposition 3.1. Suppose a 7^ yand (resp. ^) satisfies one of the following conditions: 

(1) 0(v) (resp. ^(m) = )for 6 e [0,^] (resp. $ e [0,^]), 

(2) @ (resp. ^) is Lips chitz continuous. 

Then, there exists a predictable process \]/ = (V^; j,?)i<i j<3,o<f<r in IR^^'^ such that 




for all t <T and i= 1,2. Furthermore, iWi,j,t)l<i,j<3 invertible for all t. 
Proof. Consider the process % :— e^"'{Vt + a) for t <T. Then, 

dVt = £ e-'^'OierndBi,,. 

i= 1,2,3 

(We let 01 = 02 = 01 = 0.) In other words, V is a local martingale. We first show that V is in fact a 
martingale. Suppose is Lipschitz continuous. By the Burkholder-Davis-Gundy Inequality, there 
exists a positive constant C such that 

T 

EsupVf^ < CE [ e-^"'Q{Vt)^dt 

t<T Jo 

rT rT 

< CE V^dt + C<C EV^dt + C<^ 
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by well known estimates of moments of solutions of stochastic differential equations with Lipschitz 
coefficients. On the other hand, if 0(v) = for 6 e [0, , then 

Esupy^^ < (J EVf^^ dt<C {EVtY^dt 

t<T Jo Jo 

< c£ {e'^'EVtf^ dt < c£ {e'"Vof^dt < - 

since V is a positive local martingale. Hence V is a martingale. Similarly, we can show that the 
process U, defined by Ut :— e~'^{Ut + 77) for ? < T, is a martingale when 4> satisfies Condition 1 
or 2. 

Suppose 7, a 7^ 0. Then, 

E(^£usds\^t^ = j'^Usds + E(^j^' {eWs-nds)\^)j 

t T' 

= J^Usds + J^' e^' (E {Us\^t) - n) ds 

= Usds + j^' {e^'Ut -r])ds 

= Usds + Uti J -r]{Ti-t) 

In particular. 



Gj = E {^j^' Usds\:^^ + E {^j^' Vsds\^}j - Ki 

Uo-r]Ti+( )vo-aTi-Ki 



Y J \ ^ 

rt / pYTi _ \ r-t f pUTi _ as ^ 



in which 



/e^-l\ ~ /e«^-l\ . ^ ft 

= i^-^jUo-r]Ti+i^—^jVo-aTi-Ki+ l^^J^ ¥u,sdBj,s 

/eYTi^gYtx / aTi_ at\ 

¥i,j,t = [-^) ^M^t) + [ a J 

for /= 1,2 and 7 = 1,2,3. Define Y3j,t = <yj'^tSt for j = 1,2,3. Note that \i/i,3,t = for / = 1,2. 
Since 

\ViJ,t)l<i<2,2<j<3 — y e^i-el" e"^2-e«^ j \^ Q 030(y) / 

(xT- at Y^- yt 

is invertible, so is . In the case that a (resp. 7) is equal to zero, the term ^ '^^ (resp. ^— y^) 

in the above matrix is replaced by 7} — f, and the matrix is also invertible when Ot 7^ 7. □ 
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Remark. The fact that the matrix i// can be explicitly obtained and shown to be invertible is the 
main benefit of using variance swaps to complete the market. Similar calculations for non-linear 
contingent claims like put and call options on the stock or the realized variance would have been 
much more difficult, if not impossible, to obtain. As a result, such non-linear contingent claims 
would make the hedging much more difficult in practice. Note that the processes U and V need 
not be martingales under the risk neutral measure, i.e. a, 7 7^ 0. Consequently, Jq'^^cIs is not a 
martingale and Gi t 7^ Jo '^^ds — Kj for z = 1 , 2. If that were the case, one of the two variance swaps 
would be redundant. 

From now on, we assume that « 7^ 7 and that and ^» satisfy the assumptions of the previous 
proposition. 

The next lemma implies that the best way of trading is always to use FV continuous s.f.t.s. to 
avoid liquidity costs coming from the quadratic variation of X. In this sense, trades should always 
be done at the quoted price 5(?,0). Note that even though some of the liquidity costs in Equation 
12.71 are eliminated when using continuous FV strategies, liquidity risk has not been completely 
eliminated from the model since the integral J^^ X^_dMu is still present. That is the main difference 
between our setup and the CJP model. 

If S is a special semimartingale with canonical decomposition S = N +A, i.e. in which A'^ is a 
local martingale and A is a predictable and FV process, then the norm of S is defined as 

II'5|I^2-IIvWIIl2 + II r\dAML2- 

Jo 

Lemma 3.2. Let S be a special semimartingale and X be predictable and integrable with respect 
to S. There exists a sequence of bounded continuous processes with finite variation such 

thatXQ —Xj — andX" converges to X in J^^. In particular, J X"dS / XdS in J^^. 

Proof. The statement is proved in the proof of Lemma 4. 1 of ^etin et al. ||6l. □ 

We will see that, because of the quadratic variation term in the equation of s.f.t.s., it is not 
possible to replicate exactly in general. Since continuous processes with finite variation have 
zero quadratic variation, the previous lemma will prove to be useful for the replication problem. 
Following ^etin et al. [6], we make the following definition. 

Definition 3.3. H G can be approximately replicated if there exists a sequence {X'\x"iY")n>i 
of s.f.t.s. such that Yj -^H inL^. 

In the presence of trade impacts, the process implicitly depends on X and its value at the 
maturity is 

S^t^ = St + 2X C Mu-dX,, + 2X r d[M,X]^ = ST-2X C X,,^dMu 
Jo Jo Jo 

when Xj = and Xq = 0. (Here we use the time T+ and Xt = to make sure that the hedging 

strategy is liquidated before the payoff is calculated to avoid discrepancies between the observed 

asset price before and after the maturity.) The true replication problem involves finding a s.f.t.s. 

(X, y) that replicates a terminal condition which itself depends on X. Instead, for each .x; G M, we 

consider the replication of the terminal condition given by xh{Sj^) with '■= St — 2X Jq xXu-dMu 

in which X is the solution of the replication problem in the case A = 0, £ = and x = 1. Jarrow 

[fTSl used a similar approach and interpreted Xt as the market's perception of the option's "delta" 

Xt. In the expression for Sj, x denotes the number of units to be replicated. Hence, the proposed 
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approximation for the true delta for the replication of x units is xXt. Proposition 14.21 in the next 
section gives an upper bound of the error introduced by this approximation. Let us begin by giving 
an overview of the replication problem in this simplified setting. 

3.1. Contingent Claims Replication Without Trade Impact and Liquidity Costs. When A = 
and £ = 0, the s.f.t.s. (X^., Ys)t<s<T that replicates a payoff H E satisfies 

H = Yt+ C Xu-dSu+ y r Z,u-dGi,u. (3.1) 

Also, 5° = S. 

First, note that Equation l3.1l is equivalent to the following linear backward stochastic differential 
equation (BSDE): 

3 nT 

Yt = {(yj^sXsSs + Xi,s¥ij,s + X2,s¥2j,s)dBj^„ (3.2) 

j=iJt 

0<t<T. Setting 

Zj,t = Xl,s¥l,J,s + X2,sV2,j,s+XsY3J,s (3.3) 

for j = 1 , 2, 3, the BSDE can be written as 



3 T 

Yt = H-f, I Zj,,dBj,, {0<t<T). 



(3.4) 



When H e L^, BSDE [331 has a unique solution (2,7) in x (see Pardoux and Peng (TSl 
for example). Since is invertible, we can define Xs = ^]xV ^ ^ ^) Xs by inverting 
Equation l3.3[ Then {Xs,Xs,Ys)t<s<T is the solution of l3.1[ 

3.2. The Replication Problem With Liquidity Risk. From now on, we denote by {X,x,Y) the 
solution of 13.21 with terminal condition H = h{ST)- Recall that := Sj — 2xX Jq XudMu. The 
main result of this section is the following theorem. 

Theorem 3.4. Let h : M+ ^Rbe Lipschitz continuous. Then xh{Sj) can be approximately repli- 
cated for all X eM. 

Proof. Let L > and N > 0. Let .jc G M and h satisfy the conditions of the theorem, and define 
h^iy) = h{y) if \y\ < N and h^{y) ~ h{N) otherwise. Since h is continuous on [— A^,A^], is 
bounded. Denote this bound by C/^. Define = xh^{Sj) and 

Tl = inf{0 < M < r : < or E„ > L or E,, < -^j. 
Consider the following BSDE: 

Y, = H''^'^-J\,dSs + xJ\^dM,-Y^J\udG,,, (3.5) 



for < ? < Tl in which H^'^ = E [HjlJ^n] ■ It can be re- written as 

//^.^ = Yt-Xj\l^A,du + Y,pZ,,,dB,,, 



(3.6) 
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with 

Z,-„ = a,E„S„X„-0,C(M,)x2+ £ nj,tXj,u (3.7) 

7=1,2 

for / = 2, 3, Zi u = OiLuSuXu and A„ = in which 

and ^(jc) =e4>(r"i(x))V(r"^(x)). Note that thechange of variable from (X,;t:i,;t:2) to (Zi, 22,23) 
is one-to-one because y/t is invertible. Since ^^^^^ is bounded on [0, Tl] and //^'^ G U°{j^xr), there 

exists a pair (2, y)o<r<T£ of predictable processes satisfying BSDE[3]6]by Theorem 2 of Briand 
and Hu [5J. Extend these processes to [0, T] by setting Yt = Y^^^ and 2f = for t > Tl- 

Define X and % in terms 2 with Equation 13 .71 For m > 0, define X'" = ^l{|x|<m} ^iid similarly 
for Furthermore, let 2 be given by Equation 13.71 with X and ;f replaced by X and By 
Lemma [X2I there exists a sequence {(X'"'", ^'"'")}n of bounded continuous processes with finite 
variation converging to (X'", x"") in Define 2'"'" in terms of (Y™'", X™'"). then 2'"'" ^ 2™ 
in J^^asn^ 00. Since jZ"''"dB ^ /Z^Jfi, we also have that //^ |2r'"|2j^ ^ |Zf in L\ 
Letting 

= - A f\zZfKdu + 1 rzi:dB,u and 
= 70 - A f\zlSKdu + 1 rKudB.,u, 

JO ;=l"'0 

we find 7 ;2 ^ ^-rl in ^ 

as n — > 00. Furthermore, (X ' ,x'"'"iY ' ) is a s.f.t.s. since it satisfies 
Equation [33] and = [^™'",^"^'"] = 0. 

Since 7!^ ^ //^'^ as m ^ 00^ we can find a sequence (X"'^'^,;(^"'^'^,7"'^'^)„>1 of s.f.t.s. for 

each L and such that Y"^'^ H^-^ = E (//^| in . 

Since E {Hj\^t:j^ as L — > 00 a.s. by martingale convergence, we also have convergence 

in by the Dominated Convergence Theorem. Finally since Hj converges to xh{Sj) when N 

goes to infinity, we can easily find a s.f.t.s. sequence {X'\x'\Y")n>i such that Yj — > xh{Sj) in 
as n — > 00. □ 

The economic interpretation of Theorem l3.4l is that the availability of variance swaps for trading 
makes the market approximately complete in the sense that any contingent claim with a Lipschitz 
payoff function can be approximately replicated. 

4. Analytical properties of the approximate solutions 

In the presence of price impacts, the replicating cost of x units of a contingent claim is not in 
general x times the replicating cost of 1 unit. When h be a Lipschitz continuous function, recall 
that for each x an approximating s.f.t.s. for the approximate replication of xh{Sj) is obtained from 
the solution of BSDE 13. 5[ which we denote by (X-'^.x^^Y^) to emphasize the dependence on x, 

with the terminal condition E (^xh^{S-}') for and L large. The theorems in this section give 

analytical properties of these approximate solutions for fixed L and A^. To alleviate the notation, we 
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omit the L's and A'^'s in all the expressions in this section (e.g. x = XL,h = h^, etc ...) when there is 
no possible confusion. For each t <T and each .x: G R, we define Ht{x) = ^Y/^ as the replicating cost 
per unit forx units of the claim with payoff function h. Furthermore, we let Ht{0) = lim;|.^o^f(-^)- 
The next theorem states that this limit exists and is given by the solution of the replication problem 
without trade impacts and liquidity costs of Section [BTI Recall that {X,x,Y)o<t<T denotes the 
solution of the BSDE [311 with terminal condition 

Theorem 4.1. Ht{0) =Yt = E(/z(5r) |^/) and ^X"" X in L^{dQ x dt) asx^O. 

Proof. For each x, we let (Z'^,y^)o<r<T£ be the solution of BSDE [X6l with terminal condition 
E (^xh{Sj') ■ Using the notation of the proof of Theorem 13.41 we have that A„ is bounded 
on [0, t], which means there exists a constant C > such that Ak(Zj < C|Z^p. Take < 

< Cf^ we know by the maximum principle (see 



1 



41CC^ 



First note that since 



[fT6l. Proposition 2.1) that |y/| < \x\Cn <^iora\\0<s< T. Let H'' = E (h{S^^] 
proof of Theorem 13. 41 we have shown that 



In the 



thus 



x\2 



x\H'') 



= {Y^f - 2 (^xa,{zisy: - \ \zi\^^ du + 2 j\::zidB, 

> {Yn' + / V - 2XCY,:) iZll'du + 2 Y^ZldB, 



- iYn^+l'\\z:\^du+2j\:z:dB,. 

We have that E{Jj'' A^iZf J^dul^t) 

< E{J\\Z^\^du\^t) < 2CE{x^{H''f\^t) < 2Cx^C; 



.2/-2 

N 



by taking expectations. Since 7/ ~ E (xH^ + /^^ Aj,(Z-[ J^du 



, we find 



|V/-E(//^|^,)| < x2CCjj. 



(4.1) 



Since h{S-}) ^ hiSr) a.s. asx^O, wehavethatE(//^|^f) =E(/z(5^)|^f) converges to E(/z(5r)|=^?; 
a.s. as X ^ by the Dominated Convergence Theorem. Letting x go to zero in Equation 14. 1[ we 
have //f(0) = I; = E(/j(5r)|^/). 

For the second part of the theorem, let (Z, F)o<r<r be the solution of 

3 .T 

Z = h{ST)-Y, ^j,sdBj,s {0<t<T). 

;— 1 Jt 



7=1 ■ 
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ThenZj^s={ojI-sXsSs + Xi,sWi,j,s + X2,sW2,j,s) for7= l,2andZi,, = (7iL,ZA. Moreover, E /q^ | ^Z^ - 
= E\E{h{ST)\^,)-H'\^- (Eh{ST)-^Yij\2XEj\,^{ZlJ\Y^-xt)du 



< E|/j(Sr)-/i(5|)|2 + 



4ACC 



N 



E 



X 



which goes to as x — > 0. Recall that \Xf\ = -^^^ < Zf^ on [0,t] and that the same 

inequality holds for Xt and Zi ^ for any < ? < T. Thus we find that ^X^ converges to X in L? {dQ x 
dt). ' □ 



The next proposition gives an estimate of the error introduced by using instead of S'^. 

_ 2 



Proposition 4.2. Ifh is Lipschitz continuous then E 



= 0(x^) as X ^ 0. In particular, 



E 



^0{x^)asx^O. 



Proof. In terms of Z^, the process 5*^ can be decomposed as 



5r + 



/o o'lLjS'j 
since ZJ = for s outside [0, t], whereas 



(JiLsSi 



xZigds + 



xZisdBi^s 



In the proof of the previous theorem, we found 
/o X 



El \-Zl-Zu\^du < E|/^(5r)-M'^r)l'+^^ 



E 



\Z'^.?-du 



< 2Xx^E 



fXudMu 
JO 



+ SXCC^x=0{x) 



as jc — > 0. Then, for some positive constant C, 



E 



< CE / \xZu-Z^\ du 
Jo 

< x^f{x) 



in which f{x) — 0{x) as jc — > 0. 



□ 



Under the additional assumption that h is differentiable we have that Ht{x) is also differentiable 
at X = and its derivative can be computed in terms of the solution of the replication problem 
without trade impacts. The interpretation of ///(O) is analogous to the liquidity premium per share 
Mf of the stock. It gives the additional cost per unit for the replication of the contingent claim due 
to illiquidity when the number of units replicated is small, i.e. Ht{x) ~ Ht{0) +if/(0)x when x is 
small. This is comparable to the price of the stock per share ^^(jc) = 5^(0) +MtX. 
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Proposition 4.3. Let 0<t <%. Ifh is dijferentiable everywhere except at a finite number of points, 
then Ht{x) is a.s. dijferentiable at x = and 



H[{Q) = X^(^j\{M,)X^ds\^}j-2X^(^\ST)\{s,<N}{j' tsdM, 



Proof. For x>0 small enough, we have that 



1 



+2Xe( h'{ST)l{s,<N}{I'XsdM, 



A^Zj ^ds 



< — E 

X 



2 



) +2XE(h'iST)l{Sr<N}(ST~ST) 



E 



Av {^] ds 



E 



AsZ^ds 



2X^ 
< — E 

X 



h'^iS-}) -h'^iSr) -h'{ST)l{Sr<N}(ST-ST) 



+X 



E 



A 

t \ X 



Z}]ds 



We know that the second term in the last expression goes to zero when x ^0. On the other hand, 
1 



lim^-(/.^^(5-^)-/.^^(5r; 



\im - (h^ (St +x CXdM,) -h^iSj 

h' {ST)\{Sr<N} j' ^sdMs 



since h^ is differentiable everywhere except at a finite number of points. Furthermore, note that 



<c 



X4M, 



in which C is the Lipschitz constant of h. We then get the result by the Dominated Convergence 
Theorem. □ 



5. Conclusion 

This paper extends the liquidity risk model of ^etin et al. [6J by hypothesizing the existence of 
a supply curve that evolves randomly in time and by studying the impact of trades on the supply 
curve. This leads to a new characterization of self-financing trading strategies and a sufficient con- 
dition for no arbitrage. We show the direct connection between stochastic volatility and illiquidity. 
As a result, contingent claims whose payoffs depend on the value of the asset can be approximately 
replicated with the use of variance swaps. The replicating costs of such payoffs are obtained from 
the solutions of BSDEs with quadratic growth. We show that the marginal cost and the liquidity 
premium of contingent claims can be easily computed from the solution of the replication problem 
without trade impacts. 
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